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ON Abstract 

r- T -i Burst oscillations, a phenomenon observed in a significant fraction of Type I (thermonuclear) X-ray bursts, involve 

the development of highly asymmetric brightness patches in the burning surface layers of accreting neutron stars. 

<""] Intrinsically interesting as nuclear phenomena, they are also important as probes of dense matter physics and the 

I strong gravity, high magnetic field environment of the neutron star surface. Burst oscillation frequency is also used 

O 

to measure stellar spin, and doubles the sample of rapidly rotating (above 10 Hz) accreting neutron stars with known 

jrt spins. Although the mechanism remains mysterious, burst oscillation models must take into account thermonuclear 

flame spread, nuclear processes, rapid rotation, and the dynamical role of the magnetic field. This review provides 

£>. a comprehensive summary of the observational properties of burst oscillations, an assessment of the status of the 
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theoretical models that are being developed to explain them, and an overview of how they can be used to constrain 
neutron star properties such as spin, mass and radius. 

1 Introduction 

Type I X-ray bursts are thermonuclear explosions triggered by unstable burning in the oceans of accreting neutron 
stars. The ocean is the low density fluid layer of light elements that builds up via accretion on top of the solid 



crust (see Chamel & Haensel 2008 for a discussion of the conditions under which a Coulomb plasma of ions will 
crystallize to form a solid). The basic cause of bursts, an imbalance between nuclear heating and radiative cooling 
in settling material, can be explained using simple one-dimensional calculations (for reviews see|Lewin, van Paradijs| 



& Taam 1993 and Bildsten 1998[). The past decade, however, has revealed major complexities that can no longer be 



understood within the context of such simple models ( Strohmayer & Bildsten 2006 1 



One such complexity is the development of asymmetric brightness patches, known as burst oscillations, in about 10% 



of the bursts observed with high time resolution instrumentation (Galloway et al. 20081. What drives this process 
remains unexplained, and requires us to consider flame spreading and other multidimensional effects. The highly 
atypical properties of burst oscillations from the accretion-powered pulsars also suggest an important dynamical role 



for the magnetic field. This review article provides an overview of our current understanding of burst oscillations: the 
techniques employed to detect and analyse them, their key observational properties, the status of theoretical models, 
and the ways in which they can be used to constrain stellar spin rates and the dense matter equation of state. 

1.1 Burst oscillations - a brief history 

Although a number of searches for periodic phenomena in X-ray bursts were carried out prior to the mid 1990s, and 



some tentative detections were claimed, none have stood the test of time (Lewin, van Paradijs & Taam 1993 Jongert 



& van der Klis 19961. Conclusive detection of strong periodic signals in thermonuclear bursts came only with the 
launch of the Rossi X-ray Timing Explorer (RXTE) on December 30th, 1995. 

Observations of the burster 4U 1728-34 in February 1996, only weeks after the start of science operations, led to the 



discovery of a strong 363 Hz signal in six X-ray bursts ( Strohmayer et al. 1996 1. The authors dubbed this phenomenon 



'burst oscillations'. Many of the properties that we now consider hallmarks of burst oscillations were apparent even 
in these first observations: high coherence, an upwards drift of ~ 1 Hz towards an asymptotic maximum frequency 
as the burst progressed, amplitudes of up to ~ 10 % root mean square (rms), and the apparent disappearance of the 
signal during the burst peak. 



Strohmayer et al. (1996 1 concluded that rotational modulation of a bright spot on the burning surface was the 
most likely explanation for the burst oscillations. This was based on three arguments: the expected evolutionary 
link between accreting neutron stars in Low Mass X-ray Binaries and the millisecond radio pulsars (with the former 



being progenitors of the latter via the spin recycling scenario, see Bhattacharya & van den Heuvel 1991 1; the ability 
of a localised hotspot to explain the observed amplitudes; and the high coherence of the oscillations seen in the burst 
tails. 

Between 1996 and 2002, burst oscillations were discovered in bursts from eight more sources (see Tablell]). Although 
there was still no independent confirmation of spin rate, the link between burst oscillation frequency and stellar 
rotation was nonetheless strengthened. The fact that the same frequency was seen in multiple bursts from any given 
source, and the high stability of the asymptotic frequency of the drifting burst oscillations, both pointed to rotational 



modulation of a bright spot that was near stationary in the rotating frame of the star (Strohmayer et al. 1998 Muno 
|et al. 2002[ ). The detection of burst oscillations during a superburst (a longer, more energetic burst thought to be 
due to unstable carbon burning) from 4U 1636-536 also supported this interpretation. The frequency matched that 



seen in the normal Type I bursts despite the different burst physics, suggesting an external clock (Strohmayer & 



Markwardt 2002) 



In October 2002, the 401 Hz accretion-powered millsecond X-ray pulsar SAX J1808. 4-3658 went into outburst, and 



Chakrabarty et al. (2003 I reported the first robust detection of burst oscillations from a source with an independent 
measure of the spin rate (analysis of a burst during a previous accretion episode had yielded a marginal detection, 
in't Zand et al. 2001[ ). Although the frequency of the burst oscillations drifted by several Hz in the burst rise, the 
frequency in the tail was within « 6 x 10 -3 Hz of the spin frequency. A second accretion-powered pulsar with burst 



oscillations at the spin frequency, XTE J1814-338, followed shortly thereafter (Strohmayer et al. 20031, confirming 
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the status of burst oscillation sources as nuclear-powered pulsars. 

Since this time, burst oscillations have been found in several additional sources (Table 111, including three more 
accretion-powered pulsars and two intermittent pulsars (sources that show accretion-powered pulsations only sporad- 
ically). While rotational modulation of a brightness asymmetry on the burning surface remains an integral part of 
all models, the root cause of such an asymmetry remains unresolved. 

2 Burst oscillation observations 

2.1 Principles of detection 

The standard analysis technique, when searching for periodic or quasi-periodic signals, is to use Fast Fourier Trans- 
forms to produce a power spectrum. For a comprehensive review of this topic, the reader is referred to|van der Klis| 



(19891. Here I summarize only the key points that are essential to an understanding of burst oscillation detection 
(this is of particular relevance to the discussion in Table [2] about tentative detections) and subsequent analysis of 
properties like frequency and amplitude. 

A series of X-ray photon arrival times with total duration T is binned to form a time series Xk(t), the number of 
photons ('counts') in time bin th (k = 1...N). The time resolution At — T/N is limited by the native time resolution 
of the instrument (typically 1 - 125 (is for RXTE, depending on data mode), but is also set by the desired Nyquist 
frequency /n y . This is the maximum frequency that can be studied for a given time resolution, and is given by 
/Ny = 1/(2 At). The power spectrum Pj at the Fourier frequencies Vj — j/T (j — 0,2, ...,N/2 where vn/2 — /Ny), 
using the standard Leahy normalisation (jLeahy et al. 1983]), is then given by 



P-l 

J ~ N. 



i 



N \ 2 / JV 

} J Xk COS 2-KVjtk J + I 2Z Xh Sm 27T2^£fc 

fc=l / \fc=l 



(1) 



where N 1 is the total number of photons. 

In the absence of any deterministic signal, the Poisson statistics of photon counting yield powers that are distributed 
as x 2 with two degrees of freedom (d.o.f.). The presence of a periodic signal at a given frequency will boost the power 
in the appropriate frequency bin. If we detect a large power, however, we must first evaluate its significance by 
computing the probability of obtaining such a power through noise alone. By this we mean the chances of getting 
such a high value of the power, in the absence of a periodic signal, due to the natural fluctuations in powers that 
are a by-product of photon counting statistics. We compute this using the known properties of the \ 2 distribution, 
taking into account the numbers of trials (see below). 

Many burst oscillation papers use the Z^ statistic, instead of the power spectrum. Although very similar to the 
standard power spectrum computed from a Fourier transform, it does not require that the photon arrival times be 
binned. It is defined as 
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where here n is the number of harmonics (most commonly taken to be 1), k — l,2,...n is the index used to sum 
over harmonics, iV 7 is the total number of photons, and j is an index applied to each photon. The phase 4>j f° r each 
photon is defined as 



= 2™ 

'to 



/ v(t)dt (3) 
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where v(t) is the frequency under consideration and tj is the arrival time of each photon relative to some reference 
time to- Frequency v(t) is most commonly taken to be constant (so that it can be taken outside the integral) but on 
occasion may be a function of time, for example if one wants to correct for smearing due to orbital Doppler shifts as 
the neutron star orbits the center of mass of the binary. In the absence of a deterministic (e.g. periodic) signal, Z% 
is distributed as \ 2 with 2n d.o.f.. An advantage of Z^ is that it offers an efficient way of summing harmonics. 

When searching for drifting, or quasi-periodic signals, it is common to average powers from neighbouring frequency 
bins. In searches for weak signals one can also stack, or take an average of, power spectra from many independent data 
segments (time windows) or bursts. This affects the number of degrees of freedom in the theoretical \ 2 distibution 



of noise powers, but the modified theoretical distribution is known (see |van der Klis 1989| for more details). 

To assess whether a periodic signal is indeed present, we first take the power spectrum and identify the frequency 
bins with the highest powers. One then computes the probability of obtaining such high powers through Poisson noise 
alone, using the properties of the x 2 distribution with the appropriate number of degrees of freedom. One must then 
take into account the number of trials that have been made - the more trials, the more likely one is to obtain a high 
power via statistical fluctuations alone. The number of trials is the product of the number of independent frequency 
bins, time windows, energy bands, bursts and (if appropriate) sources searched. Assessing the number of trials 
properly is complicated by the fact that many burst oscillation searches use overlapping time bins or oversampled 
frequency bins to maximise candidate signal power. Failure to account properly for this can result in the significance 
of a candidate signal being overestimated (or underestimated, if overlapping bins are treated as being independent). 
If the probability of such a high value of the power arising through Poisson noise alone is below a certain threshold 
after taking into account numbers of trials, then it is deemed significant, with the quoted significance referring to the 
chances of having obtained such a signal through Poisson noise alone. 

A complication comes from the fact that noise powers (the power distribution in the absence of a periodic signal) 
are not always distributed in accordance with the \ 2 that would be expected from constant Poisson noise. What we 
are searching for in an X-ray burst is a periodic signal superimposed on the overall deterministic rise and decay of a 
burst lightcurve. These variations give low frequency power and add sidebands to noise powers, boosting their level 



(for a nice discussion of these issues, see Fox et al. 20011. This problem is particularly acute in the rising phase of 
the burst when the lightcurve is changing rapidly. There may also be other noise components from astrophysical or 
detector processes, such as red noise associated with accretion, which may continue during the burst and contribute 
to the overall emission. 

Assessing the true distribution that the powers should take, in the absence of a periodic signal, is most easily done 
using Monte Carlo simulations. One can for example make a model of the burst lightcurve, by fitting the rise and 
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decay with exponential functions, and then use this as a basis to generate a large sample of fake lightcurves with 
Poisson counting statistics but without any periodic signals. Taking power spectra from these fake lightcurves then 
yields the true, most likely frequency-dependent, distribution of noise powers. One other advantage of the Monte 
Carlo method is that any pecularities of the data analysis process (such as the use of overlapping time windows) can 
be replicated precisely. As a rule of thumb, unless a candidate signal is particularly strong, or is repeated in multiple 
independent bursts or time windows (thereby boosting significance), Monte Carlo simulations should really be done 
to obtain an accurate assessment of significance. Alternatively, one can also try to fit the distribution of measured 
powers (excluding the frequency bins containing the candidate signals), however this is somewhat risky since it by 
definition excludes extreme values from the fit. 

One final factor that must be considered, for the accretion-powered pulsars with burst oscillations, is the likely 
continuation of the channeled accretion process during the burst. One must therefore ask whether the signal detected 
during the burst could be due to the accretion-powered pulsations. If the accretion process (accretion rate, and 
fractional amplitude of accretion-powered pulsations) remains unchanged during the burst, the addition of a large 
number of unpulsed photons during the burst would cause the measured fractional amplitude to fall (by a factor 
-^buiV-^Vacc where iVbur is the number of photons due to burst emission and iV acc the number due to accretion in 
the data segment being considered). In all cases so far, fractional amplitudes remain similar despite the fact that 
burst flux far exceeds accretion flux, indicating that the burst emission itself must also be pulsed (provided that 
the assumption of unchanged accretion holds). Whether this assumption is reasonable remains a matter of debate. 
Naively one might expect radiation pressure from a burst to hinder the accretion process. However the radiation 



may also remove angular momentum, increasing the accretion rate (Walker & Meszaros 1989 Miller & Lamb 1993 
|Ballantyne fc Everett 2005[ ). 

2.2 Measurable properties 

Once a burst oscillation signal has been identified and deemed statistically significant, one can measure several 
observational properties. Some properties can be calculated directly from the power spectrum or Z^ statistic (for the 
rest of this subsection, these two terms can be used interchangeably), whereas others require pulse profile modelling. 
The latter technique involves selecting a frequency model for the signal, v(t), and calculating a phase for each photon 
using Equation J3J, using this to assign each photon to a phase bin, and thereby building up a pulse profile (number 
of photons per phase bin). The main intrinsic properties of interest are frequency, amplitude, and phase lags. The 
first two can be measured from the power spectrum or using pulse profile modelling, but to obtain phase lags one 
must use pulse profile modelling since the power spectrum discards phase information. 

In terms of frequency, we are interested in both the absolute value and the properties of any drift. The latter is 
most commonly visualised using dynamical power spectra. These are computed by taking power spectra of short, 
usually overlapping, segments of data from the burst lightcurve, as illustrated in Figure ni The power spectra are 
then commonly plotted as contours, overlaid on the burst lightcurve. Several examples of bursts with oscillations, 
from different sources, are shown Figure [2] One can then see the frequency drift that occurs as the burst progresses. 
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Another property, related to frequency, that is often given is the coherence Q of the burst oscillation train. This is 
given by 



(4) 



Aj/fwhm 

where v is the frequency of the peak in the power spectrum and Az/fwhm its full width at half-maximum, obtained 
by fitting a Gaussian or Lorentzian function to the broad peak in the power spectrum. Note that Q values are often 
quoted for peaks where the frequency drift has been incorporated into the model (as in Equationpll, thereby rendering 



the peak in the power spectrum narrower, and hence more coherent (see for example Strohmayer & Markwardt 1999 1. 
In this case the frequency v that is used in computing Q is normally the asymptotic maximum frequency (Sections 
[T] and [1X3 I. 

The amplitude of the pulsations can be computed directly from the power spectrum, or from a folded pulse profile 
(see for example Muno, Ozel & Chakrabarty 2002 and Watts, Strohmayer & Markwardt 20051. When using the 



power spectrum, the root mean square (rms) fractional amplitude r is given by 



N-. 



(5) 



N-r-Nb, 
where P is the signal power, iV 7 the total number of photons, and Nt, the number of background photons (estimated 

from the period before or after the burst). Note that signal power is not the same as the measured power. Once 

a power has been deemed significant (i.e. unlikely to have occurred through noise alone), one must then determine 

the signal power that would give rise to the measured value, since measured power includes the effects of noise 



(Groth 19751. For details of how to compute signal power from measured power, see Vaughan et al. (19941, Muno, 



Ozel & Chakrabarty (2002 I or Watts, Strohmayer & Markwardt (2005 1. For sources that also have accretion-powered 



pulsations, this value is often corrected to remove the expected contribution due to continuing pulsed accretion during 
the burst. Other related quantities that can be computed are the amplitudes of any higher harmonics of the burst 
oscillation frequency (often referred to as the 'harmonic content'), and the dependence of the amplitude on photon 
energy. 

There are several factors to be wary of when using values of amplitude from the literature. Different groups use 
different definitions of amplitude: for a lightcurve that varies as C + Asin2nvt, different groups may quote, for 
the signal at frequency v, full fractional amplitude (2A/C), half fractional amplitude (A/C), or root mean square 
(rms) fractional amplitude (A/Cy/2). In addition, amplitude depends on the frequency model used. Fitting apparent 
frequency drift will increase amplitude, but this must be taken into account when testing models since the two are 
not independent. The third factor, specific to the power spectrum technique, is that amplitude is suppressed when 



frequency drifts between Fourier bins (van der Klis 19891. This effect can sometimes be seen in published dynamical 



power spectra: compare for example the dynamical power spectra in Figure 1 of Muno, Ozel & Chakrabarty (2002 1 



with those in Figure 1 of Watts et al. (20091, where an interbin response function is used to reduce the effects of 



this drop (although note that this does not affect the amplitudes in Muno, Ozel & Chakrabarty (20021, which were 
computed using folded pulse profiles). 
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Figure 1 Generating a dynamical power spectrum. The burst lightcurve (solid black line) is first split into segments 
(top left). In this case, there are three overlapping time windows, Ato, Ati and Afe- The time stamps To, Ti and Ti 
refer to the midpoints of each window. Computing a power spectrum for each segment (center left), one can see that 
both power and frequency vary with time (blue = low power, purple = medium power, red = high power). This can 
be plotted in power- frequency-time space (bottom left). The resulting contours of power (purple = medium, red = 
high) are projected onto the frequency-time plane and the resulting dynamical power spectrum overplotted on the 
lightcurve (right). 

Energy-dependent phase lags, which are of use when testing models, are computed by comparing pulse profiles 
constructed from photons in different energy bands. The degree to which the folded pulse profiles are offset is the phase 
lag. This type of analysis can also be used to compare the phases of burst oscillations to accretion-powered pulsations, 
for sources that show both phenomena. Care must be taken when doing this to correct for any contamination due to 
the continued presence of accretion-powered pulsations during the burst. However the magnitude of the correction is 
simple to estimate (see |Watts, Patruno fc van der Klis 2008) . 

Finally, one can analyse the dependence of the burst oscillation characteristics on the properties of the bursts and 
the accretion state of the star. Burst properties may include the total integrated flux (fluence), peak flux, the presence 
or absence of photospheric radius expansion, duration, and rise and decay timescales. In terms of accretion state, one 
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Figure 2 Light curves (grey) and dynamical power spectra (black) showing burst oscillations for bursts from several 
sources: two persistent accretion-powered pulsars (SAX J1808. 4-3658 and XTE J1814-338), two intermittent pulsars 
(HETE J1900. 1-2455 and Aql X-l) and two non-pulsars (4U 1636-536 and 4U 1728-34), using data from RXTE. 
The dynamical power spectra use overlapping 2 s windows, with new windows starting at 0.25 s intervals. We use 
a Nyquist frequency of 2048 Hz and an interbin response function to reduce artificial drops in amplitude as the 
frequency drifts between Fourier bins. The contours show Leahy normalized powers of 10-100, increasing in steps of 
10. The dashed lines on the pulsar plots indicate the spin frequency determined from accretion-powered pulsations. 

can consider the overall accretion rate (which can be estimated in various ways) and, for the pulsars, compare the 
properties of the burst oscillations to those of the accretion-powered pulsations. Note that burst type also depends 
strongly on accretion rate, so these two factors are not unrelated. 

2.3 Observational properties 

2.3.1 Overview of burst oscillation sources Secure detections of burst oscillations have now been made in 
17 sources (Table [IJ, including 5 persistent accretion-powered pulsars and 2 intermittent accretion-powered pulsars. 
Note that all are transient accretors, and the use of the term persistent to describe the presence of accretion-powered 
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pulsations through accretion episodes - as used here - should not be confused with its use to describe sources that 
are persistently, as opposed to transiently, accreting. All of the burst oscillation sources are in Low Mass X-ray 



Binaries (Tauris & van den Heuvel 20061 with orbital periods (where known) in the range 1.4 to 19 hours. There is 
one candidate ultracompact binary in this group, 4U 1738-34, but the orbital period for this source is not known, see 
|Galloway et al. 2010b| Ideally, for a detection to be secure, the same frequency should have been seen in multiple 
bursts from the same source. For one of the sources in Table IT] this is not the case. In this case the high statistical 
significance of the detection is based upon the fact that the same frequency was seen in multiple independent time 
bins within a single burst, very close to the spin frequency observed from accretion-powered pulsations. 

Table [2] gives details of several other sources for which burst oscillation detections have been claimed. Some of the 
detections are very marginal (below the standard 3a threshold) , and most rely on a power in a single independent time 
bin in a single burst. Also given in the table is a summary of the analysis procedure used to estimate the significance 
of the claimed detection. Comparing these procedures to the ideal outlined in Section |2.1| the significance of some 
of these results has probably been over-estimated. As such, they should be considered tentative until confirmed in a 
second burst. 

2.3.2 Conditions in which burst oscillations are detected For four of the five persistent pulsars, burst 
oscillations have been detected in all observed bursts, irrespective of burst properties or accretion state (see Table 
[llfor references) - although the persistent pulsars do not typically exhibit a wide range of accretion states. For the 
fifth persistent pulsar (IGR J17498-2921), oscillations have only been detected in the two brightest bursts. However 
the upper limits on the presence of oscillations in the weaker bursts are comparable to the amplitudes detected in 



the brighter bursts (Chakraborty & Bhattacharyya 20121, so they could well be present at the same level in these 



other bursts. For the other sources (the non-pulsars and intermittent pulsars), burst oscillations are detected in only 
a subset of bursts. 



A comprehensive analysis of burst and burst oscillation properties by Galloway et al. (2008) found that the prop- 
erties of the bursts (e.g. duration, recurrence times) where oscillations are detected are not unusual compared to 
the properties of bursts where oscillations are not detected. There does however appear to be a correlation with 
accretion state, as estimated from the color-color diagram, a plot of hard against soft X-ray colors. As accretion 
rate increases, Low Mass X-ray Binaries typically move from the top left (hard) to bottom right (soft), tracing out 



a Z-shaped pattern (Hasinger & van der Klis 19891. The persistent pulsars, including those with burst oscillations, 
tend to remain in the hard (low accretion rate) state. The sources that are not persistent pulsars tend to show burst 
oscillations when the source is in the soft (high accretion rate) state. 



This leads, as previously noted by Muno et al. (2001 1 and Muno, Galloway & Chakrabarty (2004), to an apparent 
relationship between the occurrence of burst oscillations, spin frequency, burst type and the presence of photospheric 
radius expansion (PRE). Sources with burst oscillation frequency < 400 Hz tend to have short, most likely helium- 
dominated bursts. Higher frequency sources tend to have longer bursts, that probably involve mixed hydrogen/helium 
burning. In the soft state, short bursts are less likely to have PRE whilst long bursts are more likely. Since burst 
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Table 1. Sources with confirmed detections of burst oscillations 



Source 



Frequency (Hz) 



References* 



Persistent accretion-powered pulsars with burst oscillations 



SAX J1808. 4-3658 


401 


IGR J17498-2921 


401 


XTE J1814-338 


314 


IGR J17511-3057 


245 


IGR J17480-2446 


11 



Chakrabarty et al. (2003 I 



Linares et al. (20111, Chakraborty & Bhattacharyya (20121 



Strohmayer et al. (2003 I 



Altamirano et al. (2010b i 



Cavecchi et al. (2011 1 



Intermittent accretion-powered pulsars with burst oscillations 
Aql X-l 550 

HETE J1900.1-2455 377 



Zhang et al. (1998 I 
Watts et al. (2009k* 



Burst oscillation sources without detectable accretion-powered pulsations 



4U 1608-522 


620 


SAX J1750. 8-2900 


601 


GRS 1741.9-2853 


589 


4U 1636-536 


581 


X 1658-298 


567 


EXO 0748-676 


552 


KS 1731-260 


521 


4U 1728-34 


363 


4U 1702-429 


329 


IGR J17191-2821 


294 



Hartman et al. (2003 l, Galloway et al. (2008 I 



Kaaret et al. (2002 I, Galloway et al. (20081 



Strohmayer et al. (19971 



Strohmayer et al. (1998), Strohmayer & Markwardt (20021 



Wijnands, Strohmayer & Franco (2001) 



Galloway et al. (2010a i 



Smith, Morgan & Bradt (19971, Muno et al. (2000) 



Strohmayer et al. (1996) 



Markwardt, Strohmayer & Swank (1999 I 



Altamirano et al. (2010a i 



t Second (confirmation) references are given for sources where the initial discovery rested on a single burst (but see also 
note J). 

* Although burst oscillations from this source have only been detected in a single burst, they were observed in multiple 
independent time bins. We therefore consider this detection to be secure. 



oscillations occur more often in the soft state, the low frequency sources are more likely to have oscillations in bursts 
without PRE, while high frequency sources are more likely to have oscillations in bursts with PRE. The distinction 
is however not absolute and the relationship between these various factors is clearly complex. For a more in-depth 



discussion of the apparent link between burst type and rotation rate, see Galloway et al. (2008) 



Given that burst oscillations (for sources that are not persistent pulsars) seem to occur preferentially in certain 
accretion states, one can ask whether the non-detection of burst oscillations in other sources is at all surprising. 



Galloway et al. (2008 I characterized position on the Z-shaped track in the color-color diagram using a variable S z 



with low values of S z corresponding to harder states, and high values to softer states. For 87% of bursts with 
oscillations, they found S z > 1.75, cementing the link between burst oscillations and the soft, high accretion rate 
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Table 2. Sources with single burst or tentative burst oscillation detections 



Frequency (Hz) 



XTE J1739-285 



GS 1826-34 



A 1744-361 



4U 0614 + 09 



SAX J 1748. 9-2021 



MXB 1730-335 



4U 1916-053 



XB 1254-690 



EXO 0748-676 



Kaarct ct al. (2007 I 
significance was csti 
bands, frequencies ; 
signal in the 
overlapping > 



and time ■ 
data. Tin 



- Act signal in a 4s window in the tail of one burst (of five observed). The 

lg Monte Carlo simulations, taking into account number of bursts, energy 

windows searched. iGalloway et aL (20081 , however, report no significant 

result seems very sensitive to the choice of time windows (start point, 



ndependent). 



|Thompson et al. (2005} report a ~ 4ct signal after stacking power 
windows from 3 separate bursts, fitting the measured powers to obtf 
account the number of bursts, time windows and frequency bins used t> 
the number of trials is underestimated: the authors made selections 
other time windows and fine-tuning both segment length and the nui 
Bhattacharyya ct al. (2006) report a ~ 6ct signal in a 4s window in 



from this s( 
of frequenci 
from this at 
I Strohmayer 



This significance is estimated 
irehed. |Galloway et al. (2008) < 
has been observed with high tin 



nfii 



X^ with 2 d.< 
m the detectii 
alution instru: 



spectra from multiple 0.25 s til 
in the distribution and taking ir 
) create the power spectra, 
in energy, in addition to searching 
nber of segments included. 

the rise of the one burst observed 
f. , taking into account the number 
n in this burst, but no other burst 



ents 



Markwardt & Kuulkcrs (20081 report 



al 



one burst (of two observed), 
frequency and time windows 
|Kaaret et al. (20 03i report s 
cancc was estimated using tl 
bins searched. This source h. 
spin frequency 442 Hz jAlta 
number of bursts and time w. 



This signific; 

, 4.4ct oscillati. 
ie x. distribut 



is estimated 



n the 

takin 



oling tail of 



g into 



th 2 d 
been discovered to be 
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Note. — The tabic gives the signifieance quoted by the authors for the candidate detection, and explains the methods used to derive this significance. 
As explained in Scction [2.l| estimates derived using x 2 with 2 d.o.f. rather than Monte Carlo simulations tend to overestimate significance. The 
procedure for estimating number of trials varies from paper to paper: ideally, proper allowance should be made for the number of bursts, energy 
bands, frequency bins and time windows searched. The SAX J1748. 9-2021 result is particularly informative in terms of illustrating the importance of 
using the correct number of trials. Whether to give equal weight in terms of numbers of trials to every burst searched remains a matter of debate. 
The number of RXTE proportional counter units active during observations varies (with more photons maximising detection chances), and in addition 
burst oscillations tend to occur preferentially in certain accretion states (Scctionb .3.2} . 



state. Table [3] summarizes the status for prolific bursters in their data set, for which burst oscillations have not 
been found. The sample includes one intermittent pulsar, SAX J1748. 9-2021. For most of these sources the range 
of accretion states that were observed was not sufficient to enable full characterization of the color-color diagram. 
This means that the position variable S z could not be calculated for these sources, and so it is difficult to judge 
whether the non-detection of burst oscillations is unexpected. For the two sources that do have a fully-characterized 
color-color diagram (4U 1705-44 and 4U 1746-37), however, there are bursts with S z > 1.75. Given that 57% of the 
bursts with S z > 1.75 studied by |Galloway et al. (2008]) showed burst oscillations, it is perhaps rather surprising 
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that no oscillations have been found in these two sources (although the source with the most high S z bursts, 4U 
1746-37, has very unusual bursting behaviour in all respects). A more rigorous analysis of whether the non-detection 
of oscillations in the other bursters can be attributed to their being observed in harder accretion states would be 
desirable. 



In terms of when burst oscillations occur during bursts, Galloway et al. (2008 1 show that although they can be 
observed at any point during the burst (rise, peak or tail), they are most commonly detected in the tails of bursts. 
Oscillation trains tend to be interrupted during the peaks of bursts, particularly (although by no means exclusively) 



during episodes of PRE (Galloway et al. 2008 Altamirano et al. 2010b I 



Burst oscillations are clearly a common feature of many normal Type I X-ray bursts. But what about the other 
types of thermonuclear burst? A number of systems (although none of the confirmed burst oscillation sources in 
Table [n) have shown intermediate duration bursts that last several hundred seconds (in't Zand, Jonker & Markwardt 
2007| [Falanga et al. 2008 Linares et al. 2009 Kuulkers et al. 20101. At present there are two models to explain the 
occurrence of such bursts, both involving the build-up and ignition of a thick layer of helium. Either the system is 
ultracompact, so that it accretes nearly pure helium (jin't Zand et al. 2005| |Cumming et al. 2006[) or a thick layer 



of helium is built up from unstable hydrogen burning at low accretion rates (Peng, Brown & Truran 2007 Cooper 



& Narayan 2007a). Due to their rarity, very few intermediate bursts have been observed with high time resolution 



instruments. Where they have, no significant burst oscillation signal has been found (Linares et al. 2009 Kuulkers 



et al. 20101. 



Burst oscillations have however been detected in one superburst (Strohmayer & Markwardt 20021. The source 



in question, 4U 1636-536, also shows burst oscillations at the same frequency in its Type I bursts. Superbursts, 
which last several hours, are thought to be caused by unstable carbon ignition deeper within the ocean than the 



Type I burst ignition point (Cumming & Bildsten 2001 Strohmayer & Brown 2002 Kuulkers 20041. The oscillations 



were detected over about 800 s during the peak of the 4U 1636-536 superburst, but not during the long, decaying 
tail. High time resolution data has only been obtained for one other superburst, from the ultracompact source 4U 
1820-30. Unlike 4U 1636-536, this source has not shown burst oscillations in its Type I bursts. Unfortunately an 
antenna malfunction on RXTE resulted in the loss of the high time resolution data from the peak of this superburst 



(Strohmayer & Brown 2002 1. However no oscillations are detected during the superburst tail (unpublished analysis) 



2.3.3 Frequencies The detection of burst oscillations in several persistent and intermittent accretion-powered 
pulsars, sources for which the spin frequency is known, has proven conclusively that burst oscillation frequency 
is very close to the known spin frequency. How close, however, depends on the source. Table HI summarizes the 
situation for the various persistent and intermittent accretion-powered pulsars with burst oscillations. In some cases 
the frequencies agree to within 10~ 8 Hz, while for others they are separated by a few Hz. 

Most burst oscillation sources exhibit frequency drift, with the frequency rising by a few Hz over the course of 
a burst. Although we do not have sufficient counts to resolve individual cycles, the drift can be modelled and the 



resulting signals can be highly coherent, with values of Q as high as 4000 (Strohmayer & Markwardt 1999). The 
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Table 3. Prolific bursters without oscillations. 



Source 



Number of bursts S z reported? Bursts with S z > 1.75 



2E 1742.9-2929 


81 


N 




Rapid Burster 


66 


N 




GS 1826-24 


65 


N 




Cyg X-2 


55 


N 




4U 1705-44 


47 


Y 


8 


4U 1746-37 


30 


Y 


20 


4U 1323-62 


30 


N 




SAX J1747.0-2853 


23 


N 




EXO 1745-248 


22 


N 




XTE J1710-281 


18 


N 




1M 0836-425 


17 


N 




SAX J1748. 9-2021 


16 


N 




4U 1916-053+ 


14 


N 




GX 17+2 


12 


N 




4U 1735-44 


11 


N 





t There is a tentative detection of a burst oscillation frequency for this source, see Table |2J 

Note. — Most burst oscillations are seen when sources are in a particular accretion state with color-color diagram position 
variable S z > 1.75, jGalloway et al. 2008| and see the text for more detail). This table lists all sources with more than 10 
bursts in the RXTE Burst Catalogue (for bursts up to June 3 2007, [Galloway et al. 2 008) for which burst oscillations have 
not been detected. In most cases S z could not be computed since the source had not been observed over the full range of 
accretion states necessary to fully characterize the color-color diagram. 



most comprehensive study of frequency evolution in burst oscillation trains to date is that carried out by|Muno et al.| 



(2002|. In most cases the frequency drifts upwards during the burst, by 1-3 Hz, to reach an asymptotic maximum 
(although the signal cannot always be tracked through the peak of the burst since the amplitude sometimes drops 
below the detectability threshold). For the two intermittent pulsars, the asympototic maximum is 0.5-1 Hz below 
the spin frequency (the drifts for the regular pulsars are more unusual, see Table B|. The asymptotic maximum for 
a given source appears to be very stable, with the fractional dispersion in asymptotic frequencies, measured using 
bursts separated by several years, being < 10 -3 . Although orbital motion might account for some of the dispersion, it 
cannot account for all of the observed variation. In terms of coherence, the majority of the oscillation trains studied 
evolved smoothly in frequency, and hence appeared to be highly coherent. However in about 30% of cases, evolution 
was not smooth, suggesting jumps in phase or frequency, or the simultaneous presence of two signals with very similar 



frequencies (Muno et al. 2002) 



The largest frequency drifts reported in the literature are of order a few Hz. Wijnands, Strohmayer fc Franco] 



(2001 1 reported an apparent drift of 5 Hz in a burst from X 1658-298, and Galloway et al. (2001 1 reported a drift of 
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Table 4. Comparison of burst oscillation and spin frequency. 



Source 



References 



Frequency comparison 



Persistent accretion-powered pulsars 
SAX J1808. 4-3658 



Chakrabarty et al. (2003 1 



IGR J17498-2921 |Chakraborty fc Bhattacharyyaj 



(20121 



XTE J1814-338 



Strohmay e'r et al. (2003| l , 
Watts, Strohmaycr & Mark- 
|wardt (200"5] , |Watts, Patruno| 



van der Klis (2008 I 



IGR J17511-3057 



Altamirano et al. (2010b I 



IGR J17480-2446 Cavecchi et al 



(201111 



Intermittent accretion-powered pulsars 
Aql X-l |Zhang et al. (1998| l, |Muno| 

|et al. (2002) 1, |Casella et al.| 



(2008) 



HETE J1900.1-2455 Watts et al. (2009 1 



Upwards drift of a few Hz in the burst rise, starting below the 
spin frequency and overshooting. In tail, frequency exceeds 
spin by a few mHz. 

The frequency of burst oscillations in the tail (they are not 
detected significantly in the rise) appears to be relatively sta- 
ble, and within ±0.25 Hz of the spin frequency. However the 
precise offset, and limits on frequency drift, have yet to be 
quantified. 

Frequency is identical to the spin frequency within the errors 
(~ 10~ 8 Hz), with no drifts. The exception is the brightest 
burst, which shows a 0.1 Hz downwards drift in the burst rise. 

Upwards drifts of ~ 0.1 Hz in the burst rise, starting below 
the spin frequency. There appears to be overshoot in some 
bursts, but whether this is significant has yet to be quantified. 
The frequency stabilises very close to the spin frequency in 
the tail (certainly within 0.05 Hz), but the precise offset has 
also yet to be quantified. 

Frequency is identical to the spin frequency within the errors 
(~ 10~ 4 Hz), with no drifts. 

Burst oscillations typically drift upwards in frequency during 
the burst. The maximum asymptotic frequency reached by the 
burst oscillations is 0.5 Hz below the spin frequency measured 
in the one brief episode where the source showed accretion- 
powered pulsations. 

Burst oscillations drift upwards by ~ 1 Hz during the burst, 
with the maximum frequency being ~ 1 Hz below the spin 
frequency. 



3.6 Hz in a burst from 4U 1916-053. In both cases, however, the burst oscillation signal drops below the detection 
threshold in the middle of the burst. Since frequency cannot be tracked continuously throughout the burst, it is not 
clear whether we are really seeing a single signal drifting, or artificially connecting separate signals or perhaps even 
noise peaks (since with many bursts a few outliers due to noise would not be unexpected). A rigorous analysis of the 
entire sample of burst oscillation observations, to determine limits on the size of the drifts that are compatible with 
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the overall data set, has yet to be carried out. 

In the one case where oscillations were observed in a superburst, the frequency drift was much smaller (0.04 Hz 



over 800 s), and entirely consistent with the expected Doppler shifts due to orbital motion (Strohmayer & Markwardt 



20021. The superburst oscillation frequency inferred from this measurement was ~ 0.4 Hz higher than any of the 



aysmptotic frequencies measured during regular Type I X-ray bursts from this source (see also Giles et al. 2002 1. If 
this source shows the same behaviour as that seen in the intermittent pulsars (Table Rl) it is likely that the superburst 
oscillation frequency is closer to the true spin frequency of the star than the frequency of the normal burst oscillations. 

2.3.4 Amplitudes The average amplitudes of burst oscillations (computed over the whole burst) are highly 



variable, even for single sources, but are mostly in the range 2-20% rms (Galloway et al. 20081. Higher average 



amplitudes have been measured for some bursts where the signal in the rise dominates (Strohmayer, Zhang & Swank 



1997|Strohmayer et al. 1998 Bhattacharyya & Strohmayer 2006 Galloway et al. 2008), although the error bars are 



larger since fewer photons can be accumulated in the short burst rise. For the oscillations detected during the 
superburst from 4U 1636-536, the amplitudes are ~ 1 % rms, lower than that measured for oscillations detected in 



Type I bursts from this source (Strohmayer & Markwardt 20021. 



Amplitudes can also vary substantially during a burst. The tendency for signals to disappear in the peaks of bursts 
has already been mentioned. This occurs most often, but not exclusively, during periods of photospheric radius 



expansion (Galloway et al. 2008), and is seen in both pulsars and non-pulsars. Significant variations in amplitude 



during bursts have been found for a number of sources ( Muno, Ozel & Chakrabarty 2002 I but there are exceptions 



where the data are consistent with a constant amplitude model (Watts, Strohmayer & Markwardt 2005). 

For the sources that are accretion-powered pulsars, one can compare the amplitudes of the burst oscillations to those 



of the accretion-powered pulsations (Chakrabarty et al. 2003 Strohmayer et al. 2003 Watts et al. 2009 Altamirano 



|et al. 2010b] |Cavecchi et al. 201 1| |Papitto et al. 2011| |Chakraborty fc Bhattacharyya 2012[ ). The burst oscillation 
amplitudes for the persistent pulsars are very different. However for each individual source the burst oscillation 
amplitude is within a few percent of (and mostly lower than) that the accretion-powered pulsations at the time of the 
burst. The difference between the amplitudes of the two types of pulsations (for at least some bursts) is statistically 
significant (Wat ts, Strohmayer fc Markwardt 2005[ ). 

Harmonic content (a significant signal at the first overtone of the burst oscillation frequency) has been found 
in most of the burst oscillation trains from the persistent accretion-powered pulsars, although not necessarily at 
the same level as that seen in the accretion-powered pulsations ([Chakrabarty et al. 2003| |Strohmayer et al. 2003| 



Watts, Strohmayer & Markwardt 2005 Altamirano et al. 2010b Cavecchi et al. 2011 1. Burst oscillation only sources 



(Muno, Ozel & Chakrabarty 20021, and burst oscillations from the intermittent pulsars (Muno, Ozel & Chakrabarty 



2002 Watts et al. 20091, tend not to have significant harmonic content, although see Bhattacharyya & Strohmayer 
|(2005[ ) for an exception in the rising phases of some bursts. 

Another property that has been studied in some detail is the dependence of burst oscillation amplitude on photon 
energy. For burst oscillation only sources, and burst oscillations from the intermittent pulsars, amplitude rises with 
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energy (Muno, Ozel & Chakrabarty 2003 Watts et al. 20091. For persistent pulsars where this property has been 



analysed, however, burst oscillation amplitude falls with energy (Chakrabarty et al. 2003 Watts & Strohmayer 2006) 



For sources that show accretion-powered pulsations (persistent or intermittent) burst oscillation amplitudes have the 
same energy dependence as the accretion-powered pulsations, irrespective of whether this is a rise or a fall (references 



above and Casella et al. 2008 1 . 



2.3.5 Phase offsets Searches for phase lags between pulse profiles constructed using photons in different 
energy bands (energy dependent phase lags) have been carried out for several burst oscillation sources. While in a 
few sources there are marginal detections of hard lags (high energy pulse arriving later than the low energy one), 



most burst oscillation profiles show no statistically significant phase offsets ( Muno, Ozel & Chakrabarty 2003 



Watts 



& Strohmayer 2006 Watts et al. 2009). This differs markedly from the behaviour of accretion-powered pulsations, 



which have significant soft lags (Cui, Morgan & Titarchuk 1998 Gierlihski, Done & Barret 2002 Galloway et al 



2002 Kirsch et al. 2004 Galloway et al. 2005 Gierlihski fc Poutanen 2005 Papitto et al. 2010) 



For the persistent accretion-powered pulsars, one can also investigate phase lags between the burst oscillation pulse 
profile and the accretion-powered pulse profile. For XTE J1814-338, which has no detectable frequency drifts in 
its burst oscillations (for all but the final burst, which occurs at much lower accretion rates and has quite different 



properties), the two sets of pulsations are completely phase- locked, with constant phase offset (Strohmayer et al 



2003 Watts, Patruno & van der Klis 2008). Indeed the burst oscillation profile is actually coincident (zero phase 



offset) with the low energy accretion-powered pulse profile. Spectral modelling indicates that the latter component 



originates from the neutron star surface ( Gierlihski, Done & Barret 2002 Gierlihski & Poutanen 2005 1. Phase-locking 



of pulsations has also been reported for IGR J175 11-3057 (Papitto et al. 20101, although frequency drifts in this source 



complicate this type of analysis (Altamirano et al. 2010b I 



3 Burst oscillation theories 

3.1 Relevant physics 

Burst oscillations are associated with the surface layers of the neutron star. As such there are several key pieces of 
physics that must be considered when developing models for the phenomenon. Before discussing specific models, it 
is worthwhile reviewing these factors. Some of these intrinsic properties remain unchanged during an X-ray burst, 
whilst others may evolve. 

3.1.1 Structure of the neutron star Figure [3] shows the stratified composition of the surface layers of a 
neutron star, with ignition depths for different burst types marked. When modelling how a burst develops, one needs 
to take into account the coupling between the various burning and non-burning ocean layers, the solid crust, and 
the photosphere. Coupling may be dynamical, chemical, or thermal. The layers will evolve and expand due to heat 



generation during the burst, particularly during episodes of strong photospheric radius expansion (Paczynski 1983 



Paczynski fc Anderson 1986 in't Zand fc Weinberg 20101 
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Photosphere 
p~l g/cm 3 



H ignition depth 
p~ 1 5 g/cm 3 



He ignition depth 
p~l0 6 g/cm 3 



C ignition depth 
p~ 1 8 g/cm 3 



Crust 
p~ 1 9 g/cm 3 



Figure 3 The outer layers of an accreting neutron star, showing where thermonuclear ignition takes place. Adapted 



from a figure in Lewin & Joss (1980 1 



3.1.2 Heat generation and transport Nuclear burning is unstable when the heat released by a themonuclear 
reaction causes an increase in reaction rate that cannot be compensated for by cooling, resulting in a thermonuclear 
runaway. Several such reactions are involved in X-ray bursts. Hydrogen can burn unstably via the cold CNO cycle: 
12 C(p, 7) 13 N(/3 + , ^) 13 C(p, 7) 14 N(p, 7) ls O(/3 + , z/) 15 N(p, a) 12 C, where the rate-controlling temperature dependent step 
is 14 N(p, 7) 15 0. Helium burns unstably via the temperature-dependent triple alpha reaction: 3 4 He — > 12 C + 7. For 



reviews of the unstable reactions, and their application to X-ray bursts, see Schwarzschild & Harm (19651, Hansen 



& van Horn (19751, Fujimoto, Hanawa & Miyaji (19811, Lewin, van Paradijs & Taam (19931, and Bildsten (19981 



Superbursts are thought to be caused by unstable burning of carbon ( Woosley & Taam 1976 Taam & Picklum 1978 



Brown & Bildsten 1998 Cumming & Bildsten 2001 Strohmayer & Brown 2002) 



There are however uncertainties inherent in modelling ignition and the progression of the thermonuclear reactions 
that affect our understanding of the heat generation process. Ignition conditions are known to depend on accretion 
rate, composition of accreted material, and heat flux from the deep crust (all quantities that are impossible to measure 



precisely). Burning, sedimentation, and mixing between bursts are also likely to play a role (Heger, Cumming & 



Woosley 2007 Peng, Brown fc Truran 2007 Piro fc Bildsten 2007 Keek, Langer fc in't Zand 2009). Despite 
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considerable theoretical progress in this area, there still exists a substantial mismatch between predicted and observed 



recurrence times for most bursters, indicating that our understanding of ignition remains poor (van Paradijs, Penninx 



& Lewin 1988 Cornelisse et al. 2003 Cumming 2005 Strohmayer & Bildsten 2006 Galloway et al. 2008 J, although 



there are a few exceptions ( Galloway & Cumming 2006 Heger et al. 2007 1 The precise details of the reaction chains 



are also not exactly known, despite excellent experimental work aimed at pinning this down (Schatz 20111. These 



uncertainties will affect not only ignition conditions but also heat generation and compositional changes throughout 
the burst (see for example [Fisker et al. 2007| |Cooper, Steiner fc Brown 2009||Davids et al. 2011| and |Schatz 2011) . 
Heat transport is also important in models of burst oscillation development. One-dimensional simulations show 



that radiative heat transport will be important in all bursts, with convection also playing an important role ( Woosley 



et al. 2004 Weinberg, Bildsten & Schatz 20061. Heat transport will also be important in determining how the ther- 



monuclear flame spreads around the star. Multi-dimensional calculations and simulations have shown that convection, 



turbulence, conduction, and advection may all play a role (Fryxell & Woosley 1982 Spitkovsky, Levin & Ushomirsky 



2002 Malone et al. 2011) 



3.1.3 Rotation and flows If burst oscillation frequency is a good measure of spin frequency (as it appears to 
be, see Section [2. 3. 3| , then most burst oscillation sources rotate rapidly. Is the rotation sufficient, however, to affect 
the evolution of the burst? The importance of rotation in burst dynamics can be estimated by calculating the Rossby 
number R . This is the ratio of inertial to Coriolis force terms in the Navier-Stokes equations, the hydrodynamical 



equations governing the motions of the fluid layers where the burst takes place (see for example Pedlosky 1987 1 



Ro 



U_ 



(6) 



where U is the velocity of the motion, L a characteristic length, and / = Attv s cos9 the Coriolis parameter {y s being 
the spin frequency of the star, and 9 the co- latitude). Rotation starts to have important dynamical effects once 
R < 1 and the Coriolis force becomes the dominant force. This means that rotational effects must be taken into 
consideration once lengthscales exceed U/f. As rotation rate increases, the effects are felt on shorter lengthscales. 
We can now estimate whether rotation is likely to affect burst oscillation mechanisms, given that oscillations have 
been found in stars that rotate in the range 11-620 Hz (Table III. 

Let us first consider the effect on flame spread. One of the effects of the Coriolis force is to balance pressure 



gradients that develop as the hot fluid expands, hence slowing spreading (Spitkovsky, Levin & Ushomirsky 2002 1. The 



appropriate speed for a flow driven by such pressure gradients is U — \/gH where g is the gravitational acceleration 
and H the scale height of the hot fluid. As rotation rate increases, the tendency of the Coriolis force to confine the 
spreading flame operates over ever shorter lengthscales. Significant confinement (slowing of flame spread) will occur 



once these are less than the size of the star, R, (Cavecchi et al. 2011 1. This occurs for rotation rates 



v a > 25 



HZ G 



H 



10 km 



TO 14 cms- 2 / V 10m / V R 
Clearly most of the burst oscillation sources are in the regime where this will be relevant. 



(7) 
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Rapid rotation will also affect oscillation modes (global wave patterns) that might develop in the ocean and 
atmosphere, excited by the thermonuclear burst. As above, simple estimates can be used to determine whether 
rotation will have an effect. Shallow water gravity waves (driven by buoyancy), for example, have a timescale r = 1/N, 



where N = \fg]~H is the Brunt- Vaisala, frequency (see for example Pringle & King 2007). The characteristic speed of 



such waves is U — H/r — \JgH. This means that modes with a wavelength A comparable to or larger than ^/gH / f 
will be sensitive to the effects of rotation. Large-scale global modes (those with the largest wavelengths, A ~ R) are 
therefore likely to be affected by rotation for most of the burst oscillation sources (compare to Equation Vty. 

Rotation can also affect ignition conditions. Rotation lowers the effective gravity g e s at the equator as compared 
to the poles, with the difference being given by 



Sg c g 

ffeff 



5% (500^) 



- / R 

10 km 



1.4 M @ 

M 



(8) 



where Qt is the Keplerian angular velocity at the surface of the neutron star. The reduced gravity leads to a higher 
rate of fuel accumulation at the equator: because of this ignition is likely to occur preferentially at equatorial latitudes 
for all but a small range of accretion rates (Spitkovs ky, Levin fc Ushomirsky 2002| [Cooper fc Narayan 2007b) ) . 

The expansion associated with the X-ray burst could also in principle lead to differential rotation or shearing flows. 
Calculations show that the change in thickness of the burning layers as they expand during a burst is A_R ~ 20 m 



Bildsten (19981. If angular momentum were conserved, this would lead to a spin down of order 



2v 3 AR 



AR 



10 km 



- r 2 Hz IsookJ Uns ) {^r) < fl > 

Over the course of the burst, the hot layer could in principle achieve multiple wraps of the underlying star. Zonal flows 
(latitudinal differential rotation involving east-west flow along latitude lines, as opposed to meridional flows, which 
involve north-south flow along longitude lines) may also develop due to temperature gradients as different latitudes 



ignite at different times ( Spitkovsky, Levin & Ushomirsky 2002 1 . In practice any shear flow (radial or latitudinal) 
would be attenuated by various frictional processes or shearing instabilities, for instance. However if shear flows do 
persist, they will affect the hydrodynamics of the surface layers, modifying for example the structure and stability of 



oscillation modes (Pringle & King 2007) 



3.1.4 Magnetic fields Measuring magnetic field strength for the burst oscillation sources is very difficult. For 
some of the pulsars, it has been possible to put limits on field strength by observing spin-down between accretion 
episodes and calculating the field strength necessary for this to occur via magnetic dipole radiation. Using this 
method it has been determined that the accretion-powered pulsar SAX J1808. 4-3658, for example, has a field ~ 10 8 



G (Hartman et al. 2008 Hartman et al. 2009). Most of the burst oscillation sources, however, are not pulsars. If 



we assume that this is because the magnetic field is too weak to cause channeled accretion at the observed accretion 
rates, then this sets an upper limit on field strength ~ 10 9 G (for a discussion of the techniques used in making this 



type of estimate see Psaltis & Chakrabarty 19991. It is possible, however, that the magnetic field is simply aligned 



with the rotation axis ( Chen & Ruderman 1993 1 or that pulsations are obscured ( Titarchuk, Cui & Wood 2002 



Gogiis., Alpar fc Gilfanov 20071. 
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The importance of magnetic fields on the dynamics of the burning layers can be estimated by calculating the 
ratio of the Lorentz force terms to the other force terms in the magnetohydrodynamic versions of the Navier-Stokes 
equations. The magnetic contribution is often split into two terms - a magnetic pressure and a magnetic tension (see 



for example Choudhuri 1998). ft should also be borne in mind that the field is unlikely to remain static during the 



burst. Flows set up by the burst (convection, spreading of fuel or flame front, shearing due to differential rotation) 



may lead to amplification and modification of the field ( Cumming & Bildsten 2000 Boutloukos, Miller & Lamb 2010 1 



While magnetic pressure is unlikely to be relevant compared to other pressure terms at the burning depth, it will 
dominate in the very outermost layers of the atmosphere. Magnetic tension, which can act to counterbalance pressure 



gradients associated with accumulating fuel and flame spread, may be dynamically important (Brown & Bildsten 



1998 Cavecchi et al. 20111. Magnetic fields will also introduce new types of oscillation and instability in the ocean 



layers (see Choudhuri 1998 for a general overview of magnetohydrodynamic oscillations, Heng & Spitkovsky 2009 



and Section 3.2.21. Magnetic fields will also affect conditions in the ocean layers even prior to ignition. Channeling 



of accreting material onto the magnetic poles of the neutron star may lead to a local over-density at the magnetic 
polar caps, in addition to temperature and composition gradients ( |Brown fc Bildsten 1 998). 

3.2 Current models for burst oscillations 

Current models for burst oscillations fall into two categories, hotspot models and global mode models, illustrated in 
very general terms in Figure HI As will be discussed in detail below, it seems likely that burning spreads initially 
from an ignition point, forming a spreading hotspot in the rising phase of the burst. The burning front then either 
stalls, to leave a persistent hotspot, or spreads around the star, exciting large-scale waves in the ocean. Both could 
lead to a brightness asymmetry in the ocean during the tail of the burst. We will treat each class of model in turn, 
noting that they may not be mutually exclusive. 

3.2.1 Hotspot models Hotspot models are based on the general principle that the thermonuclear burning is 
not occurring uniformly across the surface of the star, but is instead confined to a smaller region. This hotspot is 
then modulated by the star's rotation to give rise to the observed burst oscillations. Two types of hotspot model 
have been considered: spreading hotspots, generated temporarily in the rise as the burst ignites and flame spreads 
from the ignition site to engulf the star; and persistent hotspots caused by restriction of the burning to a small part 
of the stellar surface. 

In the tails of most bursts, the blackbody radius obtained from fitting the overall burst spectrum is comparable 



to the radius of the star, implying that the flame has engulfed the entire surface (see for example Galloway ct al 



20081. However it takes hours to days to accumulate fuel between bursts, and less than a second for a thermonuclear 



runaway to develop (Shara 1982 1. For ignition conditions to be met simultaneously across the stellar surface, the 
thermal state would need to be consistent to better than one part in ~ 10 4 ( Cumming & Bildsten 2000 1. The presence 



of slight asymmetries in accretion (due for example to magnetic channeling, or the presence of equatorial boundary 
layers), makes this unlikely - particularly for the accretion-powered pulsars. In the absence of a mechanism that could 
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Burning spreads from 
ignition location 



Model I: flame 
spread halts 



Model 2: 
Large-scale 
waves excited 



Figure 4 Simplified schematic of current burst oscillation models. Burst oscillations require the presence of a brightness 
asymmetry that can be modulated by the star's rotation so that the observer sees pulsations at close to the spin 
frequency. Ignition is expected to start at a point, with a flame front then spreading across the star (left panel, red 
= hot, blue = cold). This is known as a spreading hotspot. What happens then is unclear. One possibility is that 
the flame front stalls, so that only part of the ocean burns, leaving a hotspot that persists in the tail of the burst 
(top right). The other possibility is that the flame front spreads around the star and excites large scale waves (global 
modes) in the ocean (bottom right) . The peaks and troughs of the waves will differ in temperature so that the wave 
pattern gives rise to a brightness asymmetry. 



equalize the thermal state of the ocean to the required level, it is therefore assumed that ignition starts at a point 
and that a flame front then spreads across the stellar surface. 

The presence of such a spreading hotspot would provide a simple explanation for the presence of burst oscillations 
in the rising phase of bursts. Results from a small sample of bursts appear to support this picture: in a few bursts that 
have strong oscillations in the rising phase, the amplitude of oscillations falls as the blackbody radius of the overall 



burst spectrum increases (Strohmayer, Zhang & Swank 1997 Strohmayer et al. 1998 Bhattacharyya & Strohmayer 



20061. This is consistent with the idea of a flame front spreading to cover the star, since as the spot grows in size, 



the overall amplitude of pulsations will fall (Muno, Ozel & Chakrabarty 20021. There are however several problems 
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and theoretical questions that then arise. 

The first issue is that if point ignition and flame spread are occurring in all bursts, why do we not detect oscillations 
in the rising phase of all bursts? As discussed in Section ["2.3.2| oscillations are actually more common in the tails rather 
than the rising phase of bursts. One possibility is that this is a detection bias: it can be difficult to make detections 
in the rising phase of bursts given their short duration. Something that has yet to be done is a comprehensive study 
to check whether the current data, including upper limits, are consistent with the spreading hotspot model. 

One factor that will also affect detectability is the latitude at which ignition occurs. On a rapidly rotating star, 
ignition is expected near the equator for most accretion rates due to the lower effective gravity, which makes it easier 



to build up the column depth of material necessary for ignition (Spitkovsky, Levin & Ushomirsky 20021. There are 



however small ranges of accretion rate close to the boundary between stable and unstable burning where ignition is 



expected to occur at higher latitudes (Cooper & Narayan 2007b I. Whether a burst ignites at the equator or at higher 



latitudes can have a major effect on the detectability of burst oscillations (Maurer & Watts 2008 1. For sources with a 



substantial amount of magnetic channeling, ignition may also occur preferentially at the magnetic poles due to local 
overdensities or extra heating ( |Watts, Patruno fc van der Klis 2008). 

How the flame then spreads is equally important to burst oscillation detectability, since the flame must be able to 
spread in such a way that an azimuthal asymmetry can persist. The processes controlling flame spread in X-ray bursts 
have been an open question for years, with various heat transfer mechanisms including conduction, turbulence, and 



convection all thought to play a role (Fryxell & Wooslcy f982 Nozakura, Ikeuchi & Fujimoto 1984 Bildsten 1995 



Zingale et al. 20011. Spitkovsky, Levin & Ushomirsky (20021 also pointed out the importance of hydrodynamical 



effects, in particular the interaction between the Coriolis force and expansion of the hot burning material. The 
Coriolis force will act to slow down flame spread, thereby preserving a hotspot. Indeed some degree of confinement 
may be essential to get the flame front established and propagating, since rapid expansion and lateral spread of 



the hot material might otherwise cause the flame to stall (Zingale et al. 20031. Note however that Coriolis force 



confinement will not be effective for IGR J17480-2446, which has a rotation rate of only 11 Hz (see Section 3.1.3 and 
|Cavecchi et al. 2011[ ). The flame spread mechanism remains an unsolved problem with immense importance for burst 
oscillation models. 

If the flame spreads across the whole star during the rising phase, the spreading hotspot will dissipate. The fact 
that different regions of the star ignite at different times will leave some residual temperature asymmetry in the 
tail. However the magnitude of the predicted difference is too low to explain the amplitudes of the burst oscillations 



observed in burst tails (Cumming & Bildsten 20001. Some additional mechanism, or modification to the model, is still 



required. One possibility is the excitation of large-scale waves (Section 3.2.21. The other is that a hotspot survives 
because the burning front does not spread across the whole star. 

One way in which this might occur is if fuel is confined, with the prime candidate for a confinement mechanism 
being the star's magnetic field. Matter channeled onto the magnetic polar caps of the star will be prevented from 
crossing field lines (and hence flowing across the star) until the overpressure is sufficiently large to deform the field 



lines. Estimates by Brown & Bildsten (1998 1 show that for matter to remain confined in the ignition depth for 
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normal Type I bursts requires fields of at least 10 G, much higher than the fields strengths estimated for the burst 



oscillation sources (Section 3.1.41. Magnetohydrodynamical instabilities may also act to make magnetic confinement 



of fuel ineffective (Litwin, Brown & Rosner 2001 1. Magnetic fuel confinement on burst oscillation sources is therefore 



expected to be relatively unlikely. 

The alternative is that the flame front itself is confined: burning, once started, spreads some distance and then stalls. 
This requires heat transport to be impeded in some way. Several mechanisms have been explored in the literature. 
|Payne &: Melatos (2006 J) suggested that magnetic field evolution on an accreting neutron star might result in the 
development of strong magnetic belts around the equator that would then impede north-south heat transport. On 
rapidly rotating stars, a reduction in Coriolis-mediated confinement as the flame front tries to cross the equator might 



result in stalling ( Spitkovsky, Levin & Ushomirsky 2002 Zingale et al. 2003 1. Both of these mechanisms would confine 



burning to one hemisphere of the star: it is unclear whether this would lead to an asymmetry of sufficient magnitude to 



explain burst oscillation amplitudes. Recently Cavecchi et al. (2011 1 showed that the dynamical interaction between 



the expanding burning fuel and the radial magnetic field could in principle induce horizontal field components that 
might be large enough to restrict further spread. This would lead to more localized confinement. For this nrechanism 
to work the initial field can in principle be somewhat lower than for fuel confinement, ~ 10 9 G (consistent with 
estimates for two sources, XTE J1814-338 and IGR J17480-2446). If ignition starts at the magnetic pole, this would 
result in burst oscillations that are phase-locked to the accretion-powered pulsations and have only minimal frequency 



drift, consistent with observations from these two sources (Sections 2.3.3 and 2.3.51. This mechanism may not work 



however, for sources that appear to have lower magnetic fields, unless something in the burst process (such as a 



convective dynamo) acts to increase the field temporarily (Boutloukos, Miller & Lamb 20101 



A major question addressed even in the earliest papers on burst oscillations was how to explain the upwards 



frequency drift seen in most bursts (Section 2.3.31. Strohmayer, Zhang & Swank (19971 pointed out that horizontal 
motion of the spreading flame front could lead to drifts in frequency of a few Hz. However in this case rotation 
would have to set a preferred direction for spread, to explain why the frequency always tends to rise. They therefore 
suggested an alternative possibility: that of angular momentum conservation during vertical expansion of the heated 
burning layers. Rapid expansion would cause a hotspot to rotate more slowly, giving an observed frequency below 



the spin rate, with frequency rising as the layers cooled and contracted again (see Section 3.1.31. 

The expansion model depends on two key assumptions: that the burning layers can decouple from the underlying 
star; and that shearing within the burning layer itself is small enough to preserve the spot rather than causing 
it to be smeared out. The validity of these assumptions, and the magnitude of the expected frequency shift in 



first Newtonian gravity and then General Relativity, were examined by Cumming & Bildsten (2000 1 and Cumming 



|et al. (2002[ ). These authors examined a range of hydrodynamical coupling mechanisms, such as viscous or shearing 
instabilities, and found that it was in principle possible for the burning layer to remain decoupled for the timescale 
of the burst. Smearing out of the hotspot was more problematic, and keeping this coherent required convection or 
short wavelength baroclinic instabilities to operate inside the burning layer. Smearing was found to be particularly 
pronounced during episodes of photospheric radius expansion. The maximum frequency drifts predicted by these 
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studies were ~ 1 — 2 Hz. Although larger drifts have been reported, as explained in Section [2 . 3 . 3| the largest drifts in 
the literature are not clear cut, and the data set as a whole may in fact be consistent with a smaller maximum drift. 
Magnetic wind-up and shearing effects may also be important in determining the coupling and angular momentum 



transfer between the expanding layers and the underlying star (Cumming & Bildsten 2000 Lovelace, Kulkarni & 



Romanova 20071. A role for the magnetic field may explain the different drift properties seen in burst oscillations 



from the accretion-powered millisecond pulsars, where the magnetic field is known to be dynamically important. 

3.2.2 Surface modes Ignition of the burst and the subsequent spread of flame around the star may excite 
large-scale waves (global modes) in the neutron star ocean. Height differences associated with such oscillations would 
translate into hotter and cooler patches (the wave pattern), thus giving rise to differences in X-ray brightness. Non- 
axisymmetric modes could in principle be excited easily by the initial flame spread, and persist throughout the burst 
tail (depending on damping mechanisms). 

The problem of finding global modes of oscillation in oceans is one with a long and venerable history, due to 
its applications to the Earth's oceans. As such, papers on this topic use a rich and sometimes confusing mix 
of terminology when referring to different mode types. The basic hydrodynamical equations of continuity (mass 
conservation), momentum and energy, applied to the oscillations of a shallow stably stratified ocean on a rotating 



sphere, reduce to what is known as Laplace's tidal equation ( Longuet-Higgins 1968 Townsend 2005 Lou 20001. The 
restoring forces in this simple system are buoyancy (due to temperature or composition gradients) and the Coriolis 
force. 

Non-axisymmetric mode solutions to Laplace's tidal equation have dependence exp(irmp) on the azimuthal angle 
ip, where m is an integer. Given a mode with an azimuthal number m and frequency v r in the rotating frame of the 
star, the frequency v seen by an inertial observer would be 

v = mv a + v r (10) 

where v a is the spin frequency of the star, and the sign of v r is positive or negative depending on whether the mode is 
prograde (eastbound) or retrograde (westbound) respectively. The angular structure of the modes, which depend on 
spherical harmonics, are characterized by the usual integers I and m (see for example |Mathews fc Walker 19 70). The 
pattern has ra nodes in azimuth and \l — m\ nodes in latitude. For non-zero v r , the brightness pattern moves around 
the star, permitting an offset from the spin-frequency. Frequency drift can be naturally explained if v r depends on 
the thermal state of the ocean, since this will change as the burst evolves and cools through the tail. 
Laplace's tidal equation admits three families of mode solution: 

• Poincare modes: propagate in the retrograde (westbound) direction and reduce to pure gravity waves (restored 
by buoyancy alone) in the non-rotating limit. For this reason they are often referred to as g-modes in the 



neutron star ocean oscillation literature (see for example Bildsten, Ushomirsky & Cutler 1996 and Heyl 20041. 

Kelvin modes: propagate in the prograde (eastbound) direction and also reduce to pure gravity waves in the 
non-rotating limit. Unlike the Poincare modes, they are in geostrophic balance. This means that the horizontal 
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component of the Coriolis force balances the horizontal pressure gradient, see Pedlosky (19871. As a result they 



involve purely azimuthal motions, and are confined to an equatorial waveguide. 



Rossby modes: propagate in the retrograde (westbound) direction, and are also known as r-modes (Papaloizou 



& Pringle 1978 Saio 19821. In the limit of zero buoyancy, Rossby modes are driven purely by latitudinal 



variations in the Coriolis force. In the non-rotating limit they reduce to trivial zero-frequency solutions. If 
buoyancy is present then it will also act to restore these modes, and they are then often referred to as buoyant 



r-modes (see for example Heyl 2004 and Piro & Bildsten 2005b I 



Hey l (2004[ ) was the first to give serious consideration to whether such ocean modes might explain the properties 
of burst oscillations, particularly in the tails of bursts. He laid out a number of criteria that the modes would have to 
meet, including the fact that observed frequency should be a few Hz below the spin frequency of the star, and that 
low l,m modes would be needed to ensure that a highly sinusoidal and high amplitude burst oscillation would result 
(the finer the pattern, the lower the contrast across the stellar surface). On the basis of these criteria he determined 
that the I = 2, m = 1 buoyant r-mode was the most likely candidate. This retrograde mode has \v r \ of only a few Hz, 
giving an observed frequency slightly below the spin frequency. As the burning layers cool (changing the buoyancy), 
\u r \ gets smaller, so that observed frequency would tend towards the spin frequency. Poincare or g-modes were ruled 



out primarily on the grounds of their frequencies, since \v r \ is expected to be much larger than a few Hz (McDermott 



|fc Taam 1987] |Strohmayer fc Lee 1996| [BiTdsten fc Cutler 1995[ [Bildsten, Ushomirsky fc Cutler 1996] |Bildsten fc| 



Cumming 1998 Heyl 2004 Maniopoulou & Andersson 20041. Kelvin modes were ruled out on the grounds that they 



are prograde, so that frequency would fall as the ocean cooled, in contrast to the majority of observed frequency 
drifts. Both Kelvin and Poincare modes are also strongly confined to the equatorial regions on rapidly-rotating stars, 
which would reduce overall burst oscillation amplitudes (since there is no strong contrasting mode pattern across 
most of the surface). The buoyant r-mode, by contrast, occupies a wider band around the equator, giving greater 
contrast. 

In addition to providing a natural explanation for the sinusoidal shape, overall frequency, and amplitude of the burst 



oscillations, mode models predict that the amplitude of the burst oscillation should increase with photon energy (Heyl 



2005, [Lee fc Strohmayer 2005, Piro & Bildsten 2006[), in line with observations of the non-pulsars and the intermittent 



pulsars (although not those of the persistent pulsars). Heyl (20041 also argued that excitation of such a mode should 



be relatively easy, since the timescale of the mode matches the ~ls rise time of the burst associated with the flame 



spread. Narayan & Cooper (20071 and Cooper (2008) explored the issue of excitation in more detail, concluding that 
the mode would need to be unstable in order to grow to the required amplitudes in the tails of bursts. They considered 



various classical pumping mechanisms (Saio 1993) and conclude that the main driver is likely to be the e mechanism, 



which drives modes via nuclear energy generation. The role of the e mechanism in neutron star envelopes was studied 
previously by|McDermott fc Taam (1987[) and |Strohmayer fc Lee (1996) for thermally unstable envelopes. |Narayan] 



& Cooper (20071 and Cooper (20081, by contrast, consider thermally stable envelopes, appropriate for the tails of 



bursts. For He-rich explosions, instability is more likely when the burst is radiation pressure dominated. Instability 
was much less likely to develop in H-rich bursts. Suppression of the instability occurs during episodes of convection, 
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which may contribute to the observed disappearance of burst oscillations during the peaks of bright photospheric 
radius expansion bursts. 

The major problem with the buoyant r-mode model is that the frequency drift predicted as the burning layer cools 
is ~ 10 Hz, an order of magnitude larger than the observed values. At 10 9 K, a typical temperature in the peak of a 
burst, the buoyant r-mode frequency should be \u r \ > 10 Hz, drifting to \u r \ ~ 1 Hz as the layer cools to 10 8 K. Various 



solutions to this problem have been explored. Piro & Bildsten (2005a I pointed out the existence of another global 



mode of oscillation that is confined to the interface between the cool ocean below the burning layer and the underlying 
elastic crust. The observed frequency of this crust-interface mode would be ~ 4 Hz below the spin frequency of the 
star. If, as the buoyant r-mode drifted upwards in frequency, it could couple efficiently to this crustal-interface mode 
and transfer its energy via resonant conversion, this would naturally truncate the drift a few Hz short of the spin 



frequency (Piro & Bildsten 2005b). The frequency of the crustal interface mode is highly stable, so could account for 
the reported invariance in asymptotic frequencies. It could also explain why the burst oscillation frequency seen in a 
superburst is close to that seen in normal bursts, despite the differences in the burning layer, since the final frequency 
in both cases would be set by the crust-interface mode. Unfortunately detailed studies of the coupling between the 



buoyant r-mode and the crustal-interface wave indicate that resonant conversion is not effective (Berkhout & Levin 



2008). The fact that the separation between spin frequency and burst oscillation frequency is much less than 4 Hz 
for the persistent and accretion-powered pulsars (Table [4j> is also problematic for the crust-interface mode model. 
As a result, several alternative models have been explored. |Heyl (2004 1 pointed out that if the r-mode were instead 



confined to the photospheric layers, the drift would be smaller. This idea has yet to be explored in detail, but as 



pointed out by Berkhout & Levin (2008 1 magnetic forces will become important in these low density layers, and 
could have a substantial effect on mode frequencies. |Cumming (2005} explored what might happen to the buoyant 
r-modes if strong zonal flows (latitudinal differential rotation) were to develop in the ocean layers in the immediate 
aftermath of flame spread. While shear-modified modes could in principle match both frequencies and drifts, his 
study showed that high m modes were likely to grow much faster than the low m modes necessary to explain burst 
oscillation properties. Whether zonal flows can develop in the first place is also still an open question (|Spitkovsky,| 



Levin & Ushomirsky 20021. Most recently Heng & Spitkovsky (20091 have considered the effects that magnetic 



fields might have on oceanic mode structure, by solving the shallow water magnetohydrodynamical equations for an 
incompressible ocean with a purely radial magnetic field. In addition to magnetically modified versions of the normal 
families of modes described above, they find new solutions dominated by magnetic and rotational effects that they 
call magnetostrophic modes. Unlike the previously known modes, which are confined to the equatorial regions on a 
rapidly rotating star, these new modes are stronger near the poles. Both prograde and retrograde modes are possible: 
however the model predicts a downwards rather than an upwards frequency drift for the retrograde modes as the 
ocean cools. More detailed calculations, however, are required to determine whether or not magnetically modified 
modes such as these might play a role in explaining burst oscillations. 

Mode models remain an extremely promising mechanism to explain burst oscillations from the non-pulsars and the 
intermittent pulsars, although issues such as the frequency drift problem remain to be resolved. For the persistent 
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accretion-powered pulsars, they have more difficulty: existing mode models cannot explain the extremely close 
agreement between spin and burst oscillation frequency (Table B, which would require \v r \ <C 1 Hz for an m — 1 
mode; the unusual frequency drifts (or lack thereof); and the fact that burst oscillation amplitude decreases with 
photon energy. Magnetically modified modes may eventually account for these differences: alternatively sources with 
a dynamically important magnetic field may have a different burst oscillation mechanism (such as a magnetically 



confined hotspot, see Section 3.2.11 



4 Uses of burst oscillations 

Burst oscillations are an intrinsic part of the thermonuclear burst process. They can also, however, be used as tools 
to explore other areas of neutron star physics. In this Section I will expand briefly on some of these applications. 



4.1 Pulse profile modelling as a diagnostic of mass and radius 

The nature of the strong force, and the state of matter at extremes of temperature and density, is an extremely 
active field of research. Of particular interest is the transition from nucleons (neutrons and protons) to de-confined 
quarks and gluons or other more exotic states. At low densities this can be studied by experiments like the Large 
Hadron Collider (LHC) or in various Heavy Ion experiments such as the Relativistic Heavy Ion Collider (RHIC) 
and the Facility for Antiproton and Ion Research (FAIR). At higher densities, however, neutron stars are the only 



environment in the Universe where the transition can be explored ( Paerels et al. 2009 1 



Dense matter models, combined with the relativistic stellar structure equations, make predictions for mass and radii 



of neutron stars. By measuring these quantities, we can thus work back to the underlying nuclear physics (Lattimer 



& Prakash 20071. Burst oscillations are expected to be powerful tools in this regard because the emission comes from 



the surface of the neutron star, deep within the gravitational potential well. As photons propagate towards us, they 
are affected by relativistic effects that depend on mass and radius, such as gravitational redshift and light-bending. 
The pulse profile associated with the surface emission pattern is modified by these processes. By modelling the 
shape and energy-dependence of the observed pulse profile using relativistic ray-tracing algorithms, researchers seek 
to deconvolve these effects and hence extract information about mass and radius. Typically one has to fit for the 
surface emission pattern and the geometry (such as observer inclination), as well as the parameters of interest. 

The fact that stellar compactness will affect pulse profiles from a hotspot on the surface of the neutron star is 
long-established, with a series of papers exploring the effects of (for example) different metrics, rotation, and stellar 



oblateness (Pechenick, Ftaclas & Cohen 1983 Chen & Shaham 1989 Strohmayer 1992 Braje, Romani & Rauch 



2000 Cadeau, Leahy & Morsink 2005 Cadeau et al. 2007 Morsink et al. 2007). The techniques developed in these 



papers are generic to all types of hotspot and have been applied extensively to accretion-powered pulsations as well 
as burst oscillations. Following the discovery of burst oscillations, several authors started to model the lightcurves 
and pulse profiles that might be expected from both static and spreading hotspots or mode patterns on the surface 



of bursting neutron stars (Strohmayer, Zhang & Swank 1997 Miller & Lamb 1998 Weinberg, Miller & Lamb 2001 
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Muno, Qzel fc Chakrabarty 2002 Muno, Qzel fc Chakrabarty 2003 Viironen fc Poutanen 2004 Lee fc Strohmayer 



2005). Some attempts were also made at the inverse problem: that of fitting measured burst oscillation pulse profiles 



in order to obtain confidence contours on stellar compactness ( Nath, Strohmayer & Swank 2002 Bhattacharyya et al 



2005). Many more attempts have been made to obtain constraints by fitting accretion-powered pulsations, despite 



the complicating effects of the accretion funnel and shock on the emission, see Poutanen (2008) for a review. Burst 
oscillations offer in principle a somewhat cleaner system dominated by thermal emission. Although the constraints 
obtained using existing data were relatively weak, the potential of the technique to put tight bounds on mass and 



radius, given sufficient photons, was clear (Strohmayer 2004). As such pulse profile modelling of burst oscillations is 



now a major science driver for proposed future large area X-ray timing missions such as the Large Observatory for X- 



ray Timing (LOFT, Feroci et al. 2011 1 and the Advanced X-ray Timing Array (AXTAR, Ray et al. 20101. A concerted 



effort is now underway to address and resolve some of the outstanding issues that may introduce degeneracies into 
fits for mass and radius, such as the effect of changes in the surface emission pattern over the course of a burst. 



4.2 Neutron star spin 

The evidence outlined in Section [2. 3. 3| supports the identification of burst oscillation frequency as a good measure of 
neutron star spin, at least to within a few Hz. The main consequence of this has been to double the number of rapidly- 
rotating (above 10 Hz) accreting neutron stars with a known spin rate. Indeed the fastest spinning accreting neutron 
star currently known (4U 1608-522, with a spin of 620 Hz) is a burst oscillation source and not an accretion-powered 
pulsar. 

Identifying the spin distribution of the various classes of neutron star has been a longstanding research goal for 
many years, as they are an important element in our understanding of stellar and binary evolution. The discovery of 
rapidly rotating accreting neutron stars, in particular, was regarded as critical to confirming the recycling scenario 



for the formation of the millisecond radio pulsars (Alpar et al. 1982 Radhakrishnan & Srinivasan 1982 Bhattacharya 



& van den Heuvel 19911. With a larger number of stars in this class more detailed studies of population evolution 



are possible, revealing interesting discrepancies between models and observation (Tauris & van den Heuvel 2006 



Lorimer 2008 Hessels 2008 Kiziltan & Thorsett 2009 Tauris 20111. These point to problems in our understanding 



of, for example, the mass transfer process, magnetic field decay, and accretion torques. Other formation routes for 
millisecond radio pulsars, such as accretion-induced collapse of white dwarves, may also need to be invoked to explain 
the observed populations. 

The maximum spin rate that a neutron star can reach is also of great interest. A neutron star with a sub-millisecond 



spin period would place a tight and very clean constraint on the dense matter equation of state ( Lattimer & Prakash 



2007 Schaffher-Bielich 2008 Haensel et al. 20091. This is because the break-up, or mass-shedding, frequency (the 



spin rate above which a neutron star would fly apart as rotational forces overwhelm gravitational attraction) depends 
strongly on the composition of the star. This in turn depends on the poorly constrained behavior of the strong force 
at supranuclear densities. Although no current burst oscillation source (and indeed not even the most rapidly rotating 



radio pulsar, Hessels et al. 2006 1 rotates fast enough to rule out equations of state, the search for extremely rapid 
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rotators does motivate burst oscillation searches. Whether burst oscillations could ever develop, or be detectable, on 



such rapidly rotating stars remains however an open question (see Section 3.1.31. 

The current spin distribution of accreting neutron stars, as given by both the accretion-powered pulsars and 
the burst oscillation sources, hints at a maximum that is well below the spin rate suggested by simple estimates 



of accretion-induced spin-up over their lifetime ( Chakrabarty et al. 2003 Chakrabarty 2008 1 . If the evolutionary 



estimates are sound, there is a requirement for a braking mechanism to halt the spin-up. Magnetic braking (due to 
interaction between the stellar field and the accretion disk) is one possibility (|Ghosh fc Lamb 1978| |White fc Zhang| 



1997 Andersson et al. 2005). Another, which has generated a lot of excitement, is the emission of gravitational waves 



(Papaloizou & Pringle 1978 Wagoner 1984 Bildsten 19981. This has triggered major theoretical effort exploring the 



nature of possible gravitational wave generation mechanisms, such as crustal mountains ( Bildsten 1998 Ushomirsky, 



Cutler & Bildsten 2000 Melatos & Payne 2005 Haskell, Jones & Andersson 20061, or internal r-mode oscillations 



(Bildsten 1998 Andersson, Kokkotas & Stergioulas 1999 Levin 1999 Andersson et al. 2000 Andersson, Jones & 



|Kokkotas 2002] |Heyl 2"002"1 1 Wagoner 2002| [Nayyar fc Owen 2006[ |Bondarescu, Teukolsky fc Wasserman "20071 |HoT1 
| Andersson fc Haskell 20JTJ ). Gravitational wave braking could also make accreting neutron stars promising continuous 
wave sources for future gravitational wave detectors, although this is one application where knowing the spin to a very 



high degree of precision (certainly better than a few Hz) would be important (Watts et al. 20081. Uncertainty in the 
spin rate has a very negative impact on gravitational wave detection threshold, since increasing the parameter space 



increases the numbers of trials (reducing detection sensitivity, see Section 2.1 1 and increases the computational power 



needed to perform any search. Improving burst oscillation models so that we can make a more accurate diagnosis of 
spin rate would be of great benefit. 

5 Conclusions 

Burst oscillations are an intriguing part of the phenomenology of thermonuclear bursts on neutron stars, with wider 
uses in terms of measuring stellar spin and potentially constraining the dense matter equation of state. One of the 
legacies of the Rossi X-ray Timing Explorer has been a rich database of burst oscillation observations, giving us a clear 
understanding of their general properties and the conditions under which they develop. A satisfactory theoretical 
explanation of why and under what conditions burst oscillations develop, however, is still lacking. Burst evolution 
models that take full account of both the structure of the neutron star ocean and envelope, and the dynamical 
effects of magnetic fields, are needed to explore the mechanics of flame spread and the excitation of global oceanic 
modes. Such theoretical development is essential if we are take full advantage of the opportunities offered by the next 
generation of X-ray timing telescopes. 
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